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Abstract 

r^ \ Pecuhar velocities change the expansion rate of any observer moving relative to the 

r N ■ smooth Hubble flow. As a result, observers in a galaxy like our Milky Way can experience 

^*^ \ accelerated expansion within a globally decelerating universe, even when the drift velocities 

i-Q ■ are small. The effect is local, but the affected scales can be large enough to give the false 

impression that the whole cosmos has recently entered an accelerating phase. Generally, 



U^' 



O ■ peculiar velocities are also associated with dipole-like anisotropics, triggered by the fact that 

^^ \ they introduce a preferred spatial direction. This implies that observers experiencing locally 

^ ■ accelerated expansion, as a result of their own drift motion, may also find that the accelera- 

tion is maximised in one direction and minimised in the opposite. We argue that, typically, 
such a dipole anisotropy should be relatively small and the axis should probably lie fairly 
vj \ close to the one seen in the spectrum of the Cosmic Microwave Background. 
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1 Introduction 



Realistic observers do not simply follow the smooth universal expansion, but have their own 'pe- 
culiar' motions as well. The dipolar anisotropy of the Cosmic Microwave Background (CMB), in 
particular, has been interpreted as the result of our drift flow (at roughly 600 km/sec) relative 
^ . to the cosmic rest-frame. Peculiar motions are believed to be the result of structure forma- 

H \ tion and are expected to fade away as we move on to progressively larger lengths. Thus, the 

current concordance (WMAP5-normalised) cosmological model predicts drift velocities of ap- 
proximately 100 km/sec on scales larger than 50/i~^ Mpc [h is the Hubble parameter in units 
of 100 km/sec Mpc). Recent observations however seem to indicate bulk flows with significantly 
larger amplitude and scale [I]- [3], casting some doubt on the ACDM scenario ^. Here, we look 
into the potential implications of such large-scale drift motions, which are sometimes referred 
to as 'dark flows', for the kinematics of the associated observers. 

The first question is whether drifting observers in a perturbed, dust dominated Friedmann- 
Robertson- Walker (FRW) universe and those following the Hubble expansion could assign dif- 
ferent values (and signs) to their respective deceleration parameters. Whether, in particular, it 
is theoretically possible for a peculiarly moving observer to 'experience' accelerated expansion 
while the universe is actually decelerating. We find that the answer to this question is positive, 



when the pecuhar velocity field adds to the Hubble expansion. In other words, the drifting ob- 
server should reside in a region that expands faster than the background universe. Then, around 
every typical observer in that patch, there can be a section where the deceleration parameter 
takes negative values and beyond which it becomes positive again. Moreover, even small (rel- 
ative to the Hubble rate) peculiar velocities can lead to such local acceleration. The principle 
is fairly simple: two decelerated expansions (in our case the background and the peculiar) can 
combine to give an accelerating one, as long as the acceleration is 'weak' (with — 1 < q < 
- where q is the deceleration parameter) and not 'strong' (with q < —1) - see § 12.31 below. 
Overall, accelerated expansion for a drifting observer does not necessarily imply the same for 
the universe itself. Peculiar motions can locally mimic the effects of dark energy. Furthermore, 
the affected scales can be large enough to give the false impression that the whole universe has 
recently entered an accelerating phase. 

We then ask whether contemporary observations of large-scale bulk flows report drift veloc- 
ities fast enough to support the above scenario. Again, the answer is positive. Although the 
ACDM model predicts relatively weak peculiar motions, recent independent surveys seem to 
indicate larger rms bulk velocities on scales close to lOO/i^^ Mpc and beyond [I]-[3]. The anal- 
ysis of [2J, for example, finds velocities greater than 400 km/sec on the aforementioned lengths. 
Analogous results were also reached through the complementary methods of [3], while faster 
bulk flows have been reported by [1]. The latter study, in particular, suggests almost constant 
peculiar velocities of roughly 1000 km/sec extending as far out as 800 Mpc, with no sign of 
converging beyond this threshold. On smaller lengths (between 30 and 60 Mpc), on the other 
hand, the work of [5] indicated a positive variance in the local Hubble rate of up to 10%. Using 
these measurements, we flnd that drifting observers can experience locally accelerated expansion 
on scales of several hundred Mpc (perhaps even larger), depending on the density of the region 
in question. 

When all the observers in the drifting patch have peculiar velocities equal to the mean bulk 
flow of the region, they should see no preferred direction in the distribution of the deceleration pa- 
rameter within the patch. In reality, however, observers have their own peculiar motions, which 
are generally different from the average. In principle, this should trigger a dipolar anisotropy in 
the measured distribution of the deceleration parameter, analogous to the CMB dipole. Now, 
however, the rest frame is not that of the universal expansion but the one defined by the mean 
flow of the patch. In other words, there should be a preferred axis in the accelerated domain, 
with the acceleration maximised in one direction and minimised in the opposite. Nevertheless, 
when dealing with typical observers - those with peculiar velocities close to the average one, the 
aforementioned dipolar anisotropy must be small. Furthermore, the direction of the g-dipole 
should lie fairly close to that of the CMB, since they are both the result of the observers' drift 
flow. Interestingly, some recent reports seem to indicate that a 'cosmological axis', namely 
a weak dipolar anisotropy aligned close to the CMB dipole, might actually reside within the 
supernovae data [6|. 

We begin by outlining the basic features of the peculiar kinematics in linearly perturbed, 
dust-dominated FRW universes. Our first aim is to show that drifting observers measure dif- 
ferent expansion rates (and therefore different deceleration parameters) from those following 
the background Hubble fiow. We then explain how weak peculiar motions can lead to locally 
accelerated expansion within a globally decelerating universe. Section 3 provides the conditions 



for apparent universal acceleration and uses current observations to estimate the scale of the 
accelerated domain and the average value of the deceleration parameter there. In section 4 
we consider anisotropic peculiar flows and obtain expressions for the main kinematic variables. 
These formulae are then used to discuss dipole-like patterns in the distribution of the deceler- 
ation parameter, which are caused by the peculiar motion itself. Finally, we close with a brief 
summary of our results and a general discussion. 

2 Peculiar kinematics in perturbed FRW universes 

The frame of the Microwave Background, namely the coordinate system where the CMB dipole 
vanishes, is believed to coincide with the frame of the smooth Hubble expansion. This is the 
preferred cosmological coordinate system, with respect to which large-scale peculiar velocities 
can be defined and measured. 

2.1 Expansion rates 

Suppose that Ua is the reference 4- velocity of the CMB, relative to which the universe is an 
exact dust-dominated FRW model. In the perturbed spacetime, a typical observer in a galaxy 
like our Milky Way has 

Ua=Ua+Va, (1) 

with Va representing the observer's drift velocitylH Note that n^f ° = always and t;^ = Vav"^ <^ 1 
in our case. The latter constraint guarantees that 7 = (1 — v'^)^^''^ ~ 1 [7J. 

The mean kinematics of the peculiarly moving, the 'tilded', observers are determined by the 
volume scalar (0 = V""Ua) of their worldline congruence 0. Positive values for imply that the 
average separation between these observers increases and indicate expansion. In the opposite 
case we have contraction. Similarly, Q (with G = V^tta > 0) monitors the expansion of the 
universe. To first order in Va, we havqj 

G = G + i9, (2) 

where "& = D"'Va- This scalar measures the mean separation between neighbouring peculiar- flow 
lines and here it will always satisfy the 7?/0 <C 1 condition. In regions where "d is positive, the 
peculiar motion adds to background expansion and the drifting observers expand faster than 
the actual universe (i.e. Q > G). Here, we will always assume positive t?|j 

In multi-component systems each group of observers has its own time-direction. So, in our 
case, time can be measured relative to the CMB and along the tilded frame. Thus, differentiating 



^The Ua-Held is also timelike irrespective of the magnitude of the pecuUar velocity. Each frame defines its own 
time direction and 3-space. The tensors hab = dab+UaUb and hab = Qab+UaUb, with gab representing the spacetime 
metric, project orthogonal to Ua and Ua respectively. They also define the orthogonally projected covariant 
derivative operators, by means of Da = ha^Vb and Da = ha^'Vb (Va is the standard covariant derivative) 0. 

^Recall that we have assumed non-relativistic peculiar velocities, which implies that we can drop terms of 
order v'^ and higher from Eq. ((2)1 and the rest of our formulae. We also use geometrised units with c = 1 = SttG. 

^Systematically negative values of -0 are usually associated with structures that have decoupled from the 
background expansion, 'turned around' and started to collapse. Such gravitationally bound systems, however, 
have typical dimensions much smaller than those considered here. 



Eq. ([2]) with respect to time, we arrive at 

^ = e' + 19 , (3) 

where overdots and primes indicate time derivatives along the Ua and the Ua fields respectively 
(i.e. G = u°''Va@, '& = u"'V a'& and 0' = u'^VaG). Also, in what follows, we will always assume 
that 0' < and that If^/G'l <C 1. The scalars and 0' measure changes in the expansion rates 
along each of the aforementioned two time-directions and are determined by the associated Ray- 
chaudhuri equations |8|. In a perturbed dust-dominated FRW model with small drift velocities, 
the Raychaudhuri formulae in the CMB and the tilded frames are 

^' = -\®^-\p 'i^d ^ = _ie2-ip + DM,, (4) 

respectively. Note that p is the matter density relative to the Wa-field and p represents its tilded 
counterpart, with p = p to linear order in Va [9j. Also, Aa is the 4-acceleration measured by 
the drifting observer. This vector vanishes in the CMB frame (i.e. Aa = 0) but is nonzero in 
every other relatively moving coordinate system. In particular, Aa = Va + (0/3)fa) to linear 
order in Va [21 ED]- The 4-acceleration term in Eq. ^3) is central to our analysis. Its presence 
means that expressions (d^) and (jDa) are different, even when matter is pressureless dust and 
the peculiar velocities are smallo In other words, observers drifting with respect to the CMB 
frame have expansion rates different from that of the actual universe simply because of their 
relative motion. For our purposes, this fact represents the most significant theoretical deviation 
from the conventional single-fluid studies. 

2.2 Deceleration parameters 

If the expansion rate of a typical observer in a dust-dominated, almost-FRW universe differs 
from that of the Hubble flow, the deceleration parameters in the two frames should be different 
as well. One might then ask whether it is possible to experience accelerated expansion in one 
frame and decelerated in the other [11]. Whether, in particular, the peculiarly moving observer 
could measure a negative deceleration parameter, while the universe is actually decelerating. 

To investigate this possibility, we first need to write down the deceleration parameters as 
'measured' in the Ua and Ua frames. Expressed in terms of their corresponding volume scalars, 
these read 

l+e^j and q = -\l + —\, (5) 

respectively. Substituting the above into Eqs. dH), the latter recast into 

{l + q)e^ = e'^ + ^p and (1 + g)G2 = B^ + ^ ;5 - SD'^i^ . (6) 

Although these expressions already reveal that q and q are generally different, it helps to relate 
the two deceleration parameters directly. Following [IT], we recall that p = p and Aa = Va + 



*The fact that the scalars Q and Q' are generally different is also demonstrated in Eq. ^. Expressions @, 
however, provide additional information showing where and how this difference comes from. 



{Q/3)va to first order in Va- Then, we can combine Eqs. (I6|) to a single relation by eliminating 
the matter density. In particular, using the definition 1} = D'^Va, together with expression ^ 
and the (linear in Va) relation t? = D°"Va — Q^/3 (see [lOJ), we arrive at 



e) "G2I e 



l + g=(l + ^)(l + _) _ (i + _) , (7) 



where 0, 1? > and i?/0 <^ 1 always. Note that the above result can be also obtained from 
Eq. dSl) by means of ([2]) and © [llj. 

Expression ([7]) provides the deceleration parameter (g), as measured by an observer drifting 
with respect to the smooth Hubble flow of a dust-dominated FRW model, relative to that of the 
universe itself (g). Following ([7]), it is clear that q and q differ in general. Moreover, as long as the 
right-hand side of Eq. ([7|) remains below unity, positive values for q do not a priori guarantee the 
same for q. It is therefore theoretically possible for the tilded observer to experience accelerated 
expansion in a decelerating universe. This, however, can only happen locally. Indeed, assuming 
that the universe approaches an exact FRW model on large enough scales, peculiar velocities 
and their effects should die away as we move on to progressively larger lengths. In other words, 
there should always be a threshold beyond which i? = = t? and q = q. The question is how far 
out such a threshold could be (see § 13.21 later). 

2.3 Weakly and strongly accelerated expansion 

Before closing this section, it is worth noting that there are two types of accelerated expansion, 
depending on the value of the deceleration parameter. In particular, the condition — 1 < g < 
is equivalent to —@^/3 < < 0, which means that q and Q can be simultaneously negative. 
Naturally, analogous relations also hold between q, Q^ and 0' (see ([5])), as well as between 
the corresponding variables of the peculiar motion. On these grounds, it helps to distinguish 

between accelerated expansion with simply —l<q<0 and that with > (namely with 
q < —1). We will therefore associate condition —l<q<0 with weakly accelerated expansion 

and treat > (equivalently q < —1) as the necessary requirement for strongly accelerated 
expansion. 

The distinction between the aforementioned two types of acceleration is important because, 
although two decelerating expansions cannot lead to a strongly accelerating one, they can com- 
bine to give weak acceleration. This means that, although both the Hubble flow and the peculiar 
motion can be decelerating, the overall expansion of the drifting observer can be weakly acceler- 
ating. To be precise, suppose that 0' < —0^/3 and •& < — 1?^/3, to ensure that the background 
expansion and the peculiar flow are both decelerating. Then, starting from Eq. ([3]), we find that 

< —(0^ + 't?^)/3 < 0, which excludes the possibility of a strongly accelerated overall motion. 
Nevertheless, on using ([2]), the same expression recasts into 

e<-^02 + ^0^, (8) 

with 0, 1? > in our case. Clearly, this result allows for the possibility that 

-i02<0<O, (9) 



implying —l<q<0 for the tilded observer. In other words, pecuhar motions can lead to weakly 
accelerated expansion within the limits of the linear (the almost-FRW) approximation. |f| Noting 
that the supernovae results point towards weak acceleration, since they put the deceleration 
parameter close to —0.5 [13], we will focus on the —l<q<0 case for the rest of this report. 

3 Apparent acceleration due to peculiar motions 

The main theoretical principle following from our discussion so far, is that measuring a negative 
deceleration parameter in a frame drifting relative to the CMB (like that of our Milky Way for 
example) may simply be a local effect and does not necessarily imply a globally accelerating 
universe. Our next step is to examine the conditions for this to happen and to what extent 
current observations could support such a scenario. 

3.1 Conditions for local acceleration 

Observers will experience local acceleration, as a direct result of their own peculiar motion 
relative to the smooth Hubble flow, in a region where the right-hand side of Eq. ([7]) drops below 
unity. Recalling that in the absence of matter pressure q = il/2 and 3G' = — 0^[1 + {fl/2)], the 
latter expression recasts into 

where ri, O, t? > and O' < always. Then, noting that il may be seen as the effective density 
parameter of the region in question, rather than that of the universe itself, the condition for 
locally accelerated expansion reads 

\ / „q \ -2 / „Q \ 

<1. (11) 

So far, we have assumed drift velocities with u^ ^ 1 in a dust-dominated, almost-FRW 
universe. It should be noted here that the proximity to the FRW model has not followed from 
geometrical considerations imposed on the perturbed metric and its derivatives. Instead, in line 
with the covariant nature of our approach, the near-FRW claim follows from the conditions im- 
posed on the kinematic scalars and by demanding that the homogeneity scale of the universe lies 
inside the horizon. In particular, to ensure that the peculiar kinematics are always subdominant 
to the Hubble flow we have demanded that 'd/Q, \'d/Q'\ <C 1 at all times (see § 12.11 earlier). 
Although ■& is always assumed positive, however, its time derivative can (in principle at least) 
take either sign. Qualitatively speaking, positive values for '& will assist the (local) acceleration, 
while negative ones will do the opposite. We may therefore distinguish between the following 
three characteristic cases: 




^Perhaps the most direct way of demonstrating that two decelerating expansions can lead to an accelerated 
one, is by writing Eq. ((2]) as a/ a — {a/ a) + {a /a), where a, a and a are the three scale factors (with d, d > 0). 
Then, 5/3 = {'d/a) + {a/ a) + 2{d/a)(a/a), meaning that negative values for a and a do not a priori guarantee 
the same for a. Note that for simplicity we have used overdots for both time derivatives. 



Suppose that '& < and < -i?/©' <^ ■&/& (recalling that 0' < always). In this occasion, 
which one could argue that corresponds to peculiar motions with very slowly varying 
relative expansion rates, condition pT|) reduces to 

At the same time, the associated deceleration parameter (see Eq. (jlOp ) is given by 

Alternatively, we may assume that ^ < with < t?/0' ~ iD/Q <C 1. Then, the condition 
for locally accelerated expansion strengthens to 



1 ^\ /. i9'~^ 



with 

q = q^c^[l + ^nji^l + ^j -1. (15) 

Finally, when -iJ > 0, with -iJ/B' < and |-i?/e'| ~ -i?/0 < 1, condition ([II]) relaxes to 



and the corresponding value of the deceleration parameter is determined by 

The last of the aforementioned conditions is the the most favourable of the three, as far as local 
acceleration is concerned, while the second is the least favourable one. The first condition, on 
the other hand, corresponds to an intermediate situation and for this reason the coming sections 
will mainly focus on this case. 

Clearly, in regions where the right-hand side of (|12p . or of ()14p and (J16p . drops below unity, 
the drifting observers will experience locally accelerated expansion in a globally decelerating 
universe. Whether this happens or not, as well as the scale of the affected area, depends on the 
speed of the drift flow and on the density of the region in question. In particular the faster the 
peculiar motion and the lower the density, the larger the accelerated patch and the lower (more 
negative) the deceleration parameter. So, the next question is whether observations allow for 
peculiar velocities large enough to trigger such apparent acceleration. 




Figure 1: Typical observers in section A move along the ?ia-field, which drifts relative to the CMB 
frame (ua) with peculiar velocity Va- The patch A has positive i} = D^t;^ and so expands faster 
than the background universe (see Eq. ([2])). Inside region B the right-hand side of expression 
(flOl) drops below unity and there the 'tilded' observer experiences accelerated expansion [TT] . 



3.2 Locally accelerated expansion 

Consider an expanding spatial region A - see Fig. [H which largely complies with the FRW 
symmetries, but is still endowed with a weak, bulk peculiar velocity field that 'adds' to the 
background expansion (i.e. ?? > 0, with t?/0 ^ 1). Typical observers in A have peculiar velocities 
close to the bulk flow of the patch (both in magnitude and direction) and their expansion rate 
is determined by Eq. ([7|). According to the latter, the deceleration parameter measured by the 
tilded observer is generally different from that of the background universe. In particular, around 
any given observer in A there can be a region B where the right-hand side of ([7]) - see also (jlOp 
- drops below unity and where the deceleration parameter takes negative values. As we will see 
below, the size of the accelerated domain B and the value of the deceleration parameter there, 
depend on the specifics of the peculiar expansion and on the density distribution of the matter 
within section A. 

Peculiar velocities are difficult to assess directly. As a result, the speed of bulk drift flows 
is usually estimated by means of statistical methods (e.g. see [E])- The current concordance 
cosmological model, the ACDM, allows for drift flows no faster than ~ 100 km/sec on scales close 
to 100 Mpc that drop further as we move out to larger lengths. Recent independent surveys, 
however, have reported signiflcantly faster bulk velocities on substantially greater scales [I]- [3], 
putting the standard picture in doubt fl]. The surveys of p], in particular, indicate peculiar 
velocities as fast as ~ 1000 km/sec extending out to '^ 800 Mpc, and perhaps all the way to the 
Hubble radiusO These values translate into 'd/Q ~ 1/7 on scales close to 100 Mpc (assuming that 



®The surveys extend to lOO/i ^Mpc [2l[3], 300/i ^ Mpc and 500/i ^ Mpc J^, showing bulk velocities as large as 
500 km/sec and 1000 km/sec respectively. On smaller lengths (between 30 and 60 Mpc) the work of [5] suggests 



i? ~ 70 kni/secMpc today) |i| There, for an (effective) density parameter Q ~ 1/25, expressions 
dED-dlT]) lead to 

g^~-0.22, Qi ~ -0.11 and qj~-0.33, (18) 

respectively. In every case the deceleration parameter, as measured in the tilded frame, is 
negative. This means that the drifting observers experience locally accelerated expansion within 
a globally decelerating universe, as a result of their peculiar motion alone. 

As we move on to larger regions the relative magnitude of the drift velocity drops, the 
acceleration effect becomes progressively less prominent and eventually fades away. Nevertheless, 
depending on the value of the peculiar velocity and on the density of the region in question, 
one could reach as far out as several hundred (or even a few thousand) Mpc, before the smooth 
Hubble flow is restored. For instance, on lengths close to 500 Mpc, the surveys of [T] suggest 
that i?/0 ~ 1/35. There, maintaining the same value for the density parameter (i.e. Q ~ 1/25), 
Eqs. dlSD-dlTD give 

g^~ -0.036, Qt- -0-008 and g|~ -0.063, (19) 

respectivelylfl Further out, close to 1000 Mpc for example, the relative speed of the peculiar 
motion drops down to 'd/Q ~ 1/70. There, assuming that Q ~ 1/100 and using expressions 
([I3l)-([I7D, we find that 

q^ ~ -0.023 , q^ ~ -0.009 and q^ ~ -0.037 . (20) 

So, in this scenario, the accelerated expansion is only an apparent, local effect and beyond a 
certain scale-threshold the background Hubble flow is restored. To the unsuspecting observer, 
however, it may look as though the whole universe has moved into a phase of recent acceleration. 

4 Anisotropy in the g-distribution 

So far, we have focused on the average kinematics of peculiar motions, without incorporating 
anisotropics into our analysis. These can be triggered by a number of agents and in principle 
they could lead to observable anisotropy patterns in the spatial distribution of the deceleration 
parameter, as measured by the drifting observer. 

4.1 Anisotropic expansion rates 

To discuss deviations from isotropy in the spatial distribution of the deceleration parameter, 
which are themselves induced by anisotropics in the observer's peculiar motion, we first need to 



a (positive) variance in the local Hubble rate of up to 10%. With the possible exception of the last survey, there 
is currently no way of knowing whether these bulk flows are of the desired type (i.e. with i9 > 0). Nevertheless, 
in the absence of better data, we will use the above measurements to infer order-of-magnitude estimates for i9. 

^We use the approximate relation i? — Da«" — 9«"/9r" ~ 3v/r, where v is the magnitude of the mean bulk 
velocity and r the average size of the associated region. It is then straightforward to show that -d/O ~ v/Hr, 
with = 3-ff and H representing the background Hubble parameter. 

*The more modest drift velocities reported by [51[3] also lead to local acceleration, though then the accelerated 
region is generally smaller and the deceleration parameter less negative (see [TT] for numerical examples). 



consider the kinematics of anisotropically expanding media. In our case all the scales of interest 
are well within the (possible) curvature radius of the universe. Given that, we may ignore any 
general relativistic effects. Then, in the absence of vorticity, the generalised Hubble's law reads 



®a/3 — o ®^af5 + (^af5 , 



(21) 



where @ 



ap 



d(^pUa) and a, 



al3 



^y y^, ^, ^^ ^iA'^oi) ~ (0/3)(5a/3 - with Uq" = - represent the expansion 

and the shear tensors respectivelyO In addition, = d°^Ua and 5ap is the familiar Kronecker 
delta. Applied along the three principle shear eigen-directions, the above relation gives 



e„ 



e 



CFn 



(22) 



with X]a=i "^a = and X)a=i ©a = ©• ^s expected, positive shear values along a given spatial 
direction enhance the expansion, while negative ones slow it down. One can also use the diagonal 
components of the expansion tensor to define the effective scale factors in the aforementioned 
three eigen-directions by means of cia/act = Qa- Then, expression (j22p translates into 



a ^ f a\ do d-) 

- = !;[ — + — + — 
a o \ai 02 as 



(23) 



which relates the average with the direction-dependent scale factors. Furthermore, recalling that 
q = —aa/a?' is the average deceleration parameter, we define its direction dependent counterparts 
as 



la 



ai, 




(24) 



With this definition in hand, it is simple to obtain the relation between the average and 
the direction-dependent deceleration parameters. In particular, taking the time-derivative of 
Eq. (j23p . using definition (j24p and then employing some fairly straightforward algebra, we ar- 
rive at 



i + q 



(i + 9i) 



G 



+ (l+'?2) 



&2Y , ,, , _ .fQsV 



ej +(^ + ^^He^J 



(l + ..)(l + |i)%(l + ..)(l + |^)%(l + .3)(l + |^)^ 



(25) 



Note that the second equality has been obtained from the first after using relation (|22p . Also, 
in the case of zero shear anisotropy, namely when ai = a2 = cr^ = and qi = q2 = Qs = Q, the 
above relation reduces to a trivial identity. 



^Greek indices run from 1 to 3, while the Latin ones (used in the previous sections) take their values between 
and 3. Also, round brackets denote symmetrisation and square ones indicate antisymmetrisation. 
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4.2 Shear-like anisotropy 



Anisotropies in the observed spatial distribution of the deceleration parameter can appear for a 
number of reasons and in a variety of ways. For instance, the symmetry of domain A (see Fig. [1]) 
and the observer's position in it can induce anisotropy in the observed g-distribution. This type 
of anisotropy, however, is less of an issue when the patch A is considerably larger than B. In 
that case only observers 'living' near the edge of region A will be affected. Further deviation 
from isotropy can occur because the peculiar flow itself is not shear-free, or because the drifting 
observers have their own individual peculiar velocities, which are generally different from each 
other and also from the mean bulk motion of the patch. 

When shear anisotropy is included, the expansion tensor (see expression (|2ip ) associated 
with the overall motion of the drifting (the tilded) observer reads 



e 



1 



a/^ 



QSaB + ^0 



1 



(G + ■!?) 6al3 + O-a/3 



(26) 



where Gq^ = di^pu^), © = d'^Ua and -dajs 
expansion tensor of the drift motion, with i? 






{'&/'i)Siai3 + (Tap- The latter is the 
representing the associated volume scalar 



and Gap = 9(13^13) — {'&/3)6a/3 the peculiar shear (with a" a = 0). Also, as in § 14.11 before, we 
have assumed a non-rotating peculiar velocity field and set d^f^Va] to zero. Clearly, in the absence 
of anisotropy, aa/s = and the above relation reduces to Eq. ([2]). Moreover, written along the 
three principle shear eigen-directions, expression (I26p recasts into 



Gn 



1 



G + ??a = ^ (G + 1?) + C7„ = ^ © + a, 



(27) 



with J2a=i "^a = and J2a=i '^a = 3t?. We see again that positive/negative shear values along 
a given direction will enhance/reduce the overall expansion there. Taking the time derivative of 
the above, using definitions ([5^) and (j24p . together with the background relations q = Q/2 and 
G' = (GV3)(1 + n/2), we arrive at 



1 + ga 



i + \n 



^-l-t 



i + ^fi + 3^ 

G' V ^ 



(28) 



where we remind the reader that 0, may be seen as the effective density parameter of section 
A. This formula provides the deceleration parameter measured by the tilded observer, in the 
case of a generally anisotropic peculiar motion, along the principal shear eigen-directions. When 
combined with Eq. (j25p . the above leads to the average expression ([7|), thus verifying the math- 
ematical consistency of our analysis. 

As mentioned at the beginning of this section, anisotropies in the g-distribution can emerge 
from a variety of reasons. In what follows, we will use relation ()28|) to estimate apparent 
dipole-like patterns in the spatial distribution of the deceleration parameter that result from the 
observer's peculiar motion alone. 

4.3 Dipole-like anisotropy 

Of all the possible forms of anisotropy, the one most typically associated with peculiar motions is 
probably of the dipole type, caused by the fact that the drift flow introduces a preferred direction 
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to the observer's three dimensional space. This in turn gives rise to a characteristic (apparent) 
dipolar axis, relative to a given 'rest frame'. The CMB dipole, in particular, is believed to 
be such a Doppler-like effect, reflecting the peculiar motion of our Local Group relative to the 
rest-frame of the smooth Hubble flow. Similarly, differences between the peculiar velocity of an 
observer lying somewhere within section A and the mean flow of that patch, should lead to an 
analogous dipole-like pattern in the spatial distribution of the deceleration parameter in that 
region. 

Consider a typical observer within A, drifting with peculiar velocity close to the bulk flow 
of the patchl^ Without loss of generality, we may also assume that the observer moves along 
the first of the three shear eigen-directions (see ^ 14.21 before). To estimate the apparent dipolar 
anisotropy that the peculiar motion induces into the g-distribution, we will assign a positive 
value to the apparent shear {a\ ) along the observer's motion. Then, in the opposite direction, 

the corresponding apparent shear will he a\ = —a\ . Given that our observers are typical 
and that the aforementioned peculiar shear simply reflects their motion with respect to the mean 
bulk flow of the region (i.e. to the rest-frame of patch A), the associated anisotropy should be 
small. In other words, a\ /d <^ 1. Here, mainly for illustration purposes, we assume that 
a[ = -d/lB = —a[ . Then, expression (f28l) gives 

^ ~2 / a „h\ 

-1, (29) 





towards the direction of the motion and 

/ /I „Q \ -2 / A „k \ 

- 1 , (30) 

in the opposite way. Applying the above to a region of approximately 100 Mpc, where "!?/© ~ 1/7 
and ri ~ 1/25, while setting t?/G' <^ t9/0, we obtain 

gi"*"^ ~ -0.26 and g!"^~-0.18. (31) 

These results should be compared to the average value of g* ~ —0.22 obtained earlier for the 
same region and under the same assumptions (see Eq. (fTSb ) in § 13. 2p . According to ([HT|) , there 
is a small, but potentially measurable, dipole anisotropy in the spatial distribution of q along 
the direction of the peculiar flow. In particular, the universe seems to accelerate faster in one 
direction (that of the peculiar motion ~ where the apparent shear is positive). In the opposite 
direction, on the other hand, the acceleration seems to slow down by an equal amount (relative 
to the average) li^ This pattern is maintained as we move further out to larger scales, where the 



^"Throughout this report we only consider the kinematics of typical observers. By definition, these have peculiar 
velocities close, both in magnitude and in direction, to the mean bulk flow of the drifting domain. 

^^What happens is that the whole of region A drifts with respect to the Hubble expansion, while every observer 
inside A also moves relative to the mean bulk flow of that patch. The former motion is essentially the one causing 
the CMB dipole. The latter leads to a dipolar anisotropy in the g-distribution within A. For typical observers in 
A, with drift velocities close to the average one, the g-dipole should be less pronounced than its CMB counterpart. 
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effects of the peculiar motions weaken. On lengths close to 500 Mpc, for example, expressions 
(1291) and ^ give 

q[^^ ~ -0.046 and q['^ ~ -0.026 , (32) 

instead of the average g* = —0.036 (see Eq. (J19b ) in § 13.21) . Again, a small but in principle 
observable dipole anisotropy appears in the g-distributiono 

One might also wonder about the orientation of the dipolar axis. This is a considerably more 
demanding theoretical task. Nevertheless, it sounds plausible to argue that the g-dipole should 
not lie far from its CMB counterpart. This is based on the fact that both dipoles are caused by 
the observer's peculiar flow. On the other hand, we should not expect the associated dipole axes 
to coincide either, since their corresponding reference frames are different. In the CMB case 
the rest-frame is that of the smooth Hubble expansion, while here it coincides with the average 
bulk flow of section A. Recently, a number of reports seem to indicate that a small dipolar 
anisotropy, which is more-or-less aligned with the CMB dipole and is occasionally referred to as 
the "cosmological axis", might actually exist in the supernovae data [6]. 

5 Discussion 

A little more than a decade ago, observations of high-redshift supernovae indicated that our 
universe was expanding at an accelerating pace [Ijj . This conclusion was reached after applying 
the supernovae luminosity distances to the distance-redshift relation of an exact FRW model. 
The results have repeatedly returned negative values for q, indicating an accelerated expansion 
for our universe. In particular, the deceleration parameter was estimated close to -0.5. The same 
measurements also suggested that the accelerated phase was a relatively recent event, putting the 
transition from deceleration to acceleration around z = 0.5 [13]. As one looks at smaller redshifts, 
however, the picture seems to become less clear and there have been suggestions that the cosmic 
acceleration may have just peaked and started to slow down [15]. Explaining the supernovae 
results has dominated almost every aspect of contemporary cosmology. Dark energy, an unknown 
and elusive form of matter with negative gravitational energy, has so far been the most popular 
answer. There has been scepticism, however, since the dark-energy scenario lacks a satisfactory 
explanation in terms of fundamental physics and also faces an awkward coincidence problem 
(see [16j for a discussion). Possible alternatives to dark energy, include modifying General 
Relativity, introducing extra dimensions, or abandoning the Friedmann models altogether. 

The reason behind such drastic measures was that negative values for the deceleration pa- 
rameter seemed impossible in conventional FRW cosmologies, as well as in perturbed FRW 
models. This is not necessarily the case however. Peculiar motions can locally induce weakly 
accelerated expansion (with — 1 < g < - see § 12. 3p . even when the drift velocities are small 
and matter is simple pressure-free dust, namely within the limits of the linear approximation. 
Moreover, although the effect is local, the affected scales can be large enough (of the order of 
several hundred Mpc - perhaps even larger) to give the false impression that the whole universe 
has recently moved into a phase of global acceleration. 



^One can easily extend results H31[) and (|32|l to incorporate peculiar motions with li?/6'| — i?/S (see §[32 
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Technically speaking, the physics discussed here stems from the fact that the Raychaudhuri 
equations in the two coordinate systems (the CMB and that of a drifting observer) are not 
the same. The difference is in the 4-acceleration term, which vanishes in the CMB frame but 
takes nonzero values in any other relatively moving reference system. In practice, this means 
that the expansion rate of the drifting observer and that of the background Hubble flow are 
generally different. It is then theoretically possible to experience local acceleration in a globally 
decelerating universe. Alternatively, one could say that the observer's total motion is the 'sum' 
of the smooth Hubble expansion and of their own peculiar flow. Then, it is fairly easy to show 
that the two constituent motions (although both decelerating) can in principle combine to give 
weakly accelerated expansion (with —l<q<0~ see ^ I2.3p . The supernovae data seem to point 
towards weak acceleration. 

There are, of course, certain requirements that need to be fulfilled, if peculiar motions are to 
cause local acceleration. The main one is that the drift flow should add to the Hubble expansion 
(i.e. •& = D'^Va > 0). Put another way, the patch A in Fig. [1] should be expanding faster 
(even slightly faster) than the universe itself. In that case, the two decelerating expansions (the 
background and the peculiar) can locally combine to give a weakly accelerating one. One might 
then ask what is the likelihood for the Milky Way to rest within a region like A. Given that the 
scales of interest are of cosmological relevance, the chances could be fairly high (perhaps higher 
than 50%). The reason is that systematically negative values of t? are primarily associated with 
structures that have already decoupled from the background expansion, 'turned around' and 
started to collapse. Such gravitationally bound systems, however, are believed to be considerably 
smaller in size. 

The scale of the accelerated patch and the value of the deceleration parameter there, are 
decided by the speed of the drift flow and by the density of the matter in the section. In general, 
the larger the peculiar velocity and the lower the matter density, the larger the accelerated region 
and the lower (the more negative) the deceleration parameter. In our examples, we have adopted 
the peculiar velocities reported in the surveys of [I]. These suggest drift flows of approximately 
1000 km/sec extending as far out as 800 Mpc (perhaps even further out). Based on these, we 
found that the accelerated patch could reach lengths of the order of 1000 Mpc. The deceleration 
parameter was found to drop with scale, with its value ranging between —0.11 and —0.33 close 
to 100 Mpc and tending towards zero as one moves on to progressively larger lengths. As far as 
the matter density is concerned, we have assumed values in the (0.01, 0.1) range for the effective 
density parameter of the accelerated patch. The use of underdense regions, which expand faster 
than the background universe, renders certain analogies between this scenario and the large- 
void models [17j. Here, however, we do not need to reside near the centre of the void. Also, 
given that peculiar velocities are the result of structure formation, there are parallels with the 
'backreaction' scenarios as well (see |18j for a recent account). 

Drift flows are typically associated with dipolar anisotropics, triggered by the preferred 
spatial direction that the peculiar velocities introduce. The CMB dipole is believed to be an 
example of such a Doppler-like effect, caused by the motion of our Local Group relative to the 
rest-frame of the smooth Hubble expansion. Similarly, the motion of an observer in section A, 
relative to the mean bulk flow of that region, can induce a dipole-like anisotropy in the spatial 
distribution of the deceleration parameter. However, when the observer is a typical one - with 
drift velocity close to the average flow of the patch - the aforementioned dipolar anisotropy 
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should be small. With these in mind we have considered anisotropic peculiar motions, analysed 
their kinematics and used the associated equations to discuss dipole-like anisotropics. Our results 
indicate that drift flows are generally consistent with some degree of dipolar anisotropy in the 
spatial distribution of q. Moreover, if such a g-dipole exists, the axis should not lie far from 
that of its CMB counterpart. This sounds plausible, given that both dipoles are induced by the 
drift motion of the Milky Way, though relative to different rest-frames in each case. Recently, 
there has been a number of reports suggesting that a small dipolar anisotropy might actually 
reside in the supernovae data [6j. Moreover, the axis of the g-dipole appears to lie close to that 
of the CMB. In the literature, this has been largely treated as coincidental. Viewed from the 
perspective of this report, however, the alignment of the two dipolar axes seems more like a 
natural consequence rather than a mere coincidence. 
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